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A search for optical line emission from the two-photon decay of relic axions was conducted in 
the galaxy clusters Abell 2667 and 2390, using spectra from the VIMOS (Visible Multi-Object 
Spectrograph) integral field unit at the Very Large Telescope. New upper limits to the two-photon 
coupling of the axion are derived, and are at least a factor of 3 more stringent than previous upper 
limits in this mass window. The improvement follows from larger collecting area, integration time, 
and spatial resolution, as well as from improvements in signal to noise and sky subtraction made 
possible by strong-lensing mass models of these clusters. The new limits either require that the 
two-photon coupling of the axion be extremely weak or that the axion mass window between 4.5 
eV and 7.7 eV be closed. Implications for sterile-neutrino dark matter are briefly discussed also. 
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I. INTRODUCTION 

Axions are an obvious dark-matter candidate in some 
of the most conservative extensions of the standard model 
of particle physics. The magnitude of the charge-parity 
(CP) violating term in QCD is tightly constrained by 
experimental limits to the electric dipole moment of the 
neutron, presenting the strong CP problem P, [3, d, 
Fine tuning can be avoided through the Peccei-Quinn 
(PQ) mechanism, in which a new symmetry (the Peccei- 
Quinn symmetry) is introduced, along with a new pseu- 
doscalar particle, the axion. These ingredients dynami- 
cally drive the CP violating term to zero 0, H, Q . Via 
mixing with pions, axions pick up a mass, which is set 
by the PQ scale 0. 

Below a mass of 10^^ eV, axions will be produced 
through coherent oscillations of the PQ pseudoscalar, 
yielding a population of cold relics that dominate the 
dark-matter density 0, B l^- Above this mass, axions 
will be in thermal equilibrium at early times 0, S] . Un- 
less rria ^15 eV, the resulting relic density is insufficient 
to account for all the dark matter, but high enough that 
axions will be a nontrivial fraction of the dark matter Q ■ 
In either case, axions might be detectable through their 
couplings to standard-model particles. 

The couplings of the axion are set by the PQ scale 
and the specific axion model [1, 0, I^i the 

Dine-Fischler-Srednicki-Zhitnitski (DFSZ) axion model, 
standard-model fermions carry PQ charge, and so ax- 
ions couple to photon pairs both via electrically charged 
standard-model leptons and via mixing with pions 
[T2j . In hadronic axion models, such as the Kim-Shifman- 
Vainshtein-Zakharov (KSVZ) axion model [H, ax- 
ions do not couple to standard-model quarks or leptons 
at tree level. In KSVZ models, axions couple to gluons 
through triangle diagrams involving exotic fermions, to 
pions via gluons, and to photons via mixing with pions. 

Constraints to the two-hadron couplings of axions 
come from stellar evolution arguments, from the dura- 



tion of the neutrino burst from SN1987A, and from the 
upper limit to their cosmological density @, 0, [I^, [H, 
[TtI . [TsI . [Tgj . Upper limits to the two- photon coupling of 
the axion come from searches for solar axions 20], from 
upper limits to the intensity of the diffuse extragalactic 
background radiation (DEBRA) [21, 22], and from upper 
limits to x-ray and optical emission by galaxies and clus- 
ters of galaxies [2J, [23] . Recent searches for vacuum 
birefringence report evidence for the existence of a light 
boson [2I m, [23, m, m [^, though in a region of 
parameter space already constrained by null solar axion 
searches [20l. I32I. Issj . 

The two-photon coupling of the axion will lead to 
monochromatic line emission from axion decays to pho- 
ton pairs. Although the lifetime of the axion is much 
longer than a Hubble time, the dark-matter density in a 
galaxy cluster is sufficiently high that optical line emis- 
sion due to the decay of cluster axions could be detected. 
This line emission should trace the density profile of the 
galaxy cluster. Telescope searches for this emission were 
first suggested in Ref. [34]. In Ref. Q, this sugges- 
tion was extended to thermally produced axions. A tele- 
scope search for this emission was first attempted in Refs. 
[23l. [2^] . in which a null search imposed upper limits to 
the two-photon coupling of the axion in the mass window 
3 eV < ma < 8 eV. Less stringent constraints have been 
obtained in searches for decaying galactic axions [23l Is^ . 

In the past few years, high-precision cosmic microwave 
background (CMB) and large-scale-structure (LSS) mea- 
surements have become available and allowed new con- 
straints to axion parameters in this mass range. In par- 
ticular, axions in the few-eV mass range behave like hot 
dark matter and suppress small-scale structure in a man- 
ner much like neutrinos of comparable masses. Refer- 
ence |3(| shows that such arguments lead to an axion- 
mass bound ma < 1.05 eV. Still, given uncertainties 
and model dependences, it is important in cosmology to 
have several techniques as verification. For example, in 
extended, low-temperature MeV) reheating models 



2 



[3 71 . l38l . , light relics like axions and neutrinos are pro- 
duced non-thermally and have suppressed abundances, 
evading CMB/LSS bounds, but may still show up in tele- 
scope searches for axion decay lines [io^. Finally, other 
dark-matter candidates may show up in such searches; 
the sterile neutrino [111, |43, S, H] is one example, which 
we will discuss below. We are thus motivated to re- visit 
the searches of Refs. 23, 24] and see whether new tele- 
scopes, techniques, and observations may yield improve- 
ments. 

Refs. [H, preceded the advent of observations of 
gravitational lensing by galaxy clusters, however, and so 
the cluster mass density profiles assumed were not mea- 
sured directly, but derived using x-ray data and assump- 
tions about the dynamical state of the clusters. The 
constraints reported in Refs. [l^, [13] depend on these 
assumptions. Today, gravitational-lensing data can be 
used to determine cluster density profiles, independent 
of dynamical assumptions ^45,] . Thus, by using lensing 
mass maps and by applying the larger collecting areas 
of modern telescopes, cluster constraints to axions can 
both be tightened, and made robust. The high spatial 
resolution of integral field spectroscopy allows the use of 
lensing mass maps to extract the component of intraclus- 
ter emission that traces the cluster mass profile. Cluster 
mass models can be used to derive an optimal spatial 
weighting of the data, thus focusing on parts of the clus- 
ter where the highest signal is expected. 

To this end, we have conducted a search for optical 
line emission from the two-photon decays of thermally 
produced axions^. We used spectra of the galaxy clus- 
ters Abell 2667 (A2667) and Abell 2390 (A2390) obtained 
with the Visible Multi-Object Spectrograph (VIMOS) In- 
tegral Field Unit (IFU), which has the largest field of 
view of any instrument in its class [46*]. VIMOS is a spec- 
trograph mounted at the third unit (Mclipal) of the Very 
Large Telescope (VLT), part of the European Southern 
Observatory (ESO) in Chile [47j. In our analysis, we use 
mass models of the clusters derived from strong-lensing 
data, obtained with the Hubble Space Telescope (HST), 
using the Wide Field Planetary Camera #2 (WFPC-2). 

We obtain new upper limits to the two-photon coupling 
of the axion in the mass window 4.5 eV < ma < 7.7 eV 
(set by the usable wavelength range of the VIMOS IFU) 
of C < 0.003 - 0.017. The two-photon couphng of the 
axion, ^, is defined in Eq. ([3]) and discussed in Section 
im Although we search a smaller axion mass range than 
Refs. [23l[24]. our upper limits improve on past work by a 
factor of 2.1 — 7.1, depending on the candidate axion mass 
and how the limits of Ref. |23i, 2M rescaled to correct 



for today's best-fit cosmological parameters and more ac- 
curate cluster mass profiles. Our data rule out the canon- 
ical KSVZ and DFSZ models in the 4.5 eV - 7.7 eV win- 
dow. However, theoretical uncertainties in quark masses 
and pion couplings may allow for a much wider range of 
values of ^ than the canonical KSVZ and DFSZ models 
allow, as emphasized by Ref. [S], thus motivating the 
search for axions with smaller values of ^. 

A quick estimate shows that our level of improvement 
is not unexpected: The collecting area of the VLT is a 
factor of (8.1/2.1)^ greater than the 2.1m telescope at 
Kitt Peak National Observatory (KPNO) used in Ref. 
p^ . Our integration time is a factor of 10.8 ksec/3.6 ksec 
greater. The IFU allows us to cover 3.4 times as much 
of the field of view as the spectrographs used at KPNO. 
Thus we estimate that our collecting area should be a 
factor of ^ 160 higher than that of Ref. [13]. If there is 
no signal, and if we are noise limited, we would expect a 
constraint to flux that is a factor of ~ 13 more stringent 
than that of Ref. [23] , and, since oc , upper limits to 
f that are ~ 3.5 times more stringent than those reported 
in Ref. ;24]. 

We begin by reviewing the relevant theory and proceed 
to describe our observations. We then summarize our 
data analysis technique. The new limits to axion parame- 
ter space are then discussed along with other constraints. 
We conclude by pointing out the potential of conducting 
such work with higher redshift clusters. For consistency 
with the assumptions used to derive the strong-lensing 
maps used in our analysis, we assume a ACDM cosmology 
parameterized hy h — 0.71, Qni ~ 0.30, and f^A = 0.70, 
except where explicitly noted otherwise. 



II. THEORY 

To predict the expected intensity of the optical sig- 
nal due to axion decay, given the mass distribution of a 
galaxy cluster, we need to know the total mass density 
in axions. If rria > 10~^ eV, standard thermal freeze-out 
arguments show that the mass density of thermal axions 
today is 
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(1) 



where ma,eV is the axion mass in units of eV 0, H, HI] . 

Thermal axions in our mass range of interest, which 
become nonrelativistic when 4.5 eV < T < 7.7 eV, 
will have a velocity dispersion today of (^a/c^)^^^ = 
4.9 X 10~^m~gY, low enough that axions will bind to 
galaxy clusters [2^ . The maximum axionic mass frac- 
tion x"^^"" of a cluster is fUlii 



^ Based on observations made with ESO Telescopes at the Paranal 
Observatories (program ID 71.A-3010), and on observations 
made with the NASA/ESA Hubble Space Telescope, obtained 
from the data archive at the Space Telescope Institute. STScI is 
operated by the association of Universities for Research in As- 
tronomy, Inc. under the NASA contract NAS5-26555. 



(2) 



where (Tiooo is the velocity dispersion of the cluster in 
units of 1000 km s^^, and 0250 is the cluster core radius in 
units of 250 h^^ kpc. For ma.cv > 3 and typical cluster 
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values, there is ample phase space to accommodate an 
axionic mass fraction of x^^ = f2a/^^in- 

Axions decay to two photons via two channels. In one, 
axions couple to neutral pions through their two-gluon 
coupling and two QCD vertices. These pions then de- 
cay to photon pairs through QED triangle anomaly dia- 
grams. In the other mechanism, axions couple directly to 
standard-model fermions through triangle anomaly dia- 
grams, which then couple to photon pairs through QED 
vertices. The total decay rate is derived by taking the 
sum of the matrix elements for these two processes, and 
then incorporating the relevant kinematic factors, yield- 
ing an axion lifetime of [1, 

T = 6.8 X lO^'^r^m-Jv s, 

whereC= |(£^/A^-l-92±0.08). (3) 

The values of E and N depend on the axion model cho- 



sen, but by parameterizing r in terms of ^, we are able 
to obtain model independent upper limits to ^. The neg- 
ative sign comes from interference between the different 
channels for the two-photon decay of axions. The uncer- 
tainty in the theoretical value of ^ comes from uncertain- 
ties in the quark masses and pion-decay constant, and 
may in fact be larger than indicated by Eq. A com- 
plete cancellation of the axion's two-photon coupling is 
possible for KSVZ models, in which E/N ~ 2, and even 
for DFSZ axion models, in which E/N = 8/3 |4i|. It is 
thus hasty to claim that an upper limit on ^ truly rules 
out axions; it always pays to keep looking. 

The rest-frame wavelength of the axion-decay line is 
Aa = 24, 800A/ma.oVi and the line width is dominated 
by Doppler broadening. If the axion has a cosmological 
density given by Eq. ([T]) and its mass fraction in the 
cluster is Xa^ then the observer- frame specific intensity 
from axion decay is 



/a„ = 2.68 X 10^ 
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cgs. 
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where Aq denotes wavelength in the observer's rest 
frame, A^ = ^^0/(1 + ^ci), cgs denotes units of 
specific intensity (ergs cm~^ s~^ A arcsec^^), 
S'(zci) = da(2ci)/ [c/ (100 km s-i Mpc^^)] is a di- 
mensionless angular-diameter distance, and E12 = 
E/ (lO^^MQpixeP^) is the normalized surface mass den- 
sity of the cluster. If for some reason (e.g., low- 
temperature reheating [1^), the cosmological axion mass 
density is lower than indicated by Eq. ([1]), then the in- 
tensity in Eq. ^ is decreased accordingly. 

The cluster mass density was determined by fitting pa- 
rameterized potentials to the locations of gravitationally 
lensed arcs. In our mass maps, one pixel is 0.5 arcsec 
across. The intensity predicted by Eq. @ is comparable 
with that of the night-sky continuum, and so it is crucial 
to obtain a good sky subtraction when searching for an 
axion-decay line in clusters. Fortunately, the spatial de- 
pendence of the cluster density and the expected signal 
provides a natural way to separate the background from 
an axion signal, as discussed in Section IIVB I 



III. OBSERVATIONS 

A. Imaging Data 

To construct the lensing models used in our analy- 
sis and to mask out IFU fibers corresponding to clus- 
ter galaxies and other bright sources, we used images of 



A2667 and A2390 obtained with the HST and the VLT. 
The cluster A2667 was observed with HST on October 
10-11, 2001, using WFPC-2 in the F450W, F606W, and 
F814W filters, with total exposure times of 12.0 ksec, 
4.00 ksec, and 4.00 ksec, respectively [s^l. The cluster 
A2390 was observed with HST on December 10, 1994, 
using WFPC-2 in the F555W and F814W fihers and to- 
tal exposure times of 8.40 ksec and 10.5 ksec jsTi]. Af- 
ter pipeline processing, standard reduction routines were 
used with both clusters to combine the frames and re- 
move cosmic rays. Figs. [T] and [5] are images of the cluster 
cores, with iso-mass contours overlaid from our best-fit 
lensing models. 

On May 30 and June 1, 2001, near- infrared J-band 
and H-band observations of A2667 were obtained with 
ISAAC on the VLT [5jQ] . The total exposure times for the 
J and H band ISAAC data were 7.93 ksec and 6.53 ksec, 
respectively. The final seeing was 0.51" and 0.58" in the 
J and H bands, respectively. 



B. VIMOS Spectra 

The massive galaxy clusters A2667 and A2390 were 
observed with VIMOS, between June 27 and 30, 2003 
[sol l5l| . The IFU is one of three modes available on VI- 
MOS, and consists of 4 quadrants, each containing 1600 
fibers. We used an instrumental setup in which each fiber 
covered a region of 0.67" in diameter. A single pointing 
covered a 54" x 54" region of the sky. Roughly 10% of 
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FIG. 1: Image of the Abell 2667 cluster core imaged with HST 
in the F450W, F606W, and F814W filters. The white (thin 
yellow) square shows the IFU field of view, which is 54" x 54" . 
North is to the top and east is to the left. Note the strongly 
magnified gravitational arc north-east of the central galaxy. 
The white curves correspond to iso-mass contours from the 
lens model; the dark gray (red) line is the critical line at the 
redshift of the giant arc. The field of view is centered on 
aj20oo=23:52:28.4, (5j2ooo=-26:05:08. At a redshift oi z = 
0.233, the angular scale is 3.661 kpc/arcsec. 




FIG. 2: Image of the Abell 2390 cluster core imaged with 
HST in the F450W, F606W, and F814W filters. The white 
(thin yellow) squares correspond to the IFU field of view in 
different pointings. The white curves correspond to iso-mass 
contours from the lens model. The dark gray (red) line is the 
critical line at the redshift of the giant arc, labelled 1. Each 
square is 54" x 54". North is to the top and east is to the 
left. The field of view is centered on aj2ooo=21:53:36.970, 
5j2000=+17:41:44.66. At a redshift oi z = 0.228, the angular 
scale is 3.601 kpc/arcsec. 



the IFU field of view is unresponsive because of incom- 
plete fiber coverage. A low resolution blue (LR-Blue) 
grism was used, covering the wavelength range 3500A 
to TOOOA with spectral resolution R « 250 and disper- 
sion 5.355A/pixcl. The FWHM of the axion-decay line 
is 195A fJiooo "^a oV ^'^ LR-Blue grism can re- 
solve this line, without spreading a faint signal over too 
many wavelength pixels. Unfortunately, because spectra 
from contiguous pseudo-slits (sets of 400 spectra) on the 
CCD overlap, the first and last 50 pixels on most of the 
raw spectra are unusable, reducing the spectral range to 
4000A — 68OOA, corresponding to an axion mass-range 
of 4.5 < ma,cv < 7.7 at the nearly identical redshifts 
(z « 0.23) of the two clusters. Further observational 
details are discussed in Ref. [soj . 

The total exposure time for each cluster was 10.8 ksec 
(4 X 2.70 ksec exposures). Calibration frames were ob- 
tained after each of the exposures, and a spectrophoto- 
metric standard star was observed. In order to compen- 
sate for the presence of a small set of bad fibers, we used 
an offset between consecutive exposures. At a redshift 
of z = 0.233 (A2667), the IFU covers a physical region 
of 198 kpc X 198 kpc in the plane of the cluster. At a 
redshift of z = 0.228 (A2390), the IFU covers a physical 
region of 195 kpc x 195 kpc. 



C. Reduction of IFU data 

If axions exist and are present in the halos of massive 
galaxy clusters, a distinct spectral feature will appear in 
VIMOS-IFU data. At a rest-frame wavelength Aa, we 
will observe a spatially extended emission line whose in- 
tensity traces the projected dark-matter density. Reveal- 
ing such a faint, spatially extended signal requires great 
care in correcting for fiber efficiency and in subtracting 
the sky background, because the instrument itself can 
impose spatial variation in the sky background through 
varying IFU fiber efficiency. 

The VIMOS-IFU data were reduced using the VIMOS 
Interactive Pipeline Graphical Interface (VIPGI), and 
the authors' own routines [H^]. References [13, give 
both a detailed description of the methods and an as- 
sessment of the quality of VIPGI data reduction. The 
reduction steps that precede the final combination of the 
dithered exposures into a single data cube are performed 
on a quadrant by quadrant basis. The main steps are the 
following [43, [Ho, [El, 113, H^]: extract spectra from the raw 
CCD data at each pointing, calibrate wavelength, remove 
cosmic rays, determine fiber efficiencies, subtract the sky 
background, and calibrate fiux. 

The exposures were bias subtracted. Cosmic-ray hits 
were removed with an efficient automatic routine based 
on a cr-clipping algorithm, which exploits the fact that 
cosmic-ray hits show strong spatial gradients on the CCD 
[47t , in contrast to the smoother spatial behavior of gen- 
uine emission lines. In Ref. spectra were obtained 
using a limited number of long-slit exposures, so the re- 



5 



moval of a small number of incorrectly identified cosmic- 
ray hits could thwart a search for line emission from 
decaying axions. An axion-decay line, however, must 
smoothly track the density profile of the cluster. Our 
spectra are highly spatially-resolved, and so cosmic-ray 
hits can be removed safely using our cleaning algorithm. 
Using the raw CCD spectral traces, we verified that the 
signals removed by the cleaning algorithm bore the dis- 
tinctive visual signatures of cosmic-ray hits. 

VIPGI usually determines fiber efficiencies by normal- 
izing to the flux of bright sky lines; this technique yielded 
data cubes with prominent bright and dark patches (each 
covering ~ 20 x 20 fibers). It is conceivable that an acci- 
dental correlation of these patches with the cluster den- 
sity profile could lead to a spurious axion signal. To avoid 
this possibility, we measured fiber efficiency using high 
signal to noise continuum arc-lamp exposures (analogous 
to flat-fielding for images and henceforth referred to as 
flat- fielding). The resulting flat-ficlded data cubes were 
much less patchy, and were thus used for all subsequent 
analysis. 

The VIMOS IFU does not have a dedicated set of 
fibers to determine the sky-background level. VIPGI usu- 
ally determines the sky statistically at each wavelength. 
VIPGI first groups the fibers in each quadrant into three 
sets according to the shape of a user selected sky-emission 
line, and then takes the statistical mode of the counts in 
each set and subtracts it from the counts measured in 
each fiber in the set |5^. Although axions (and thus 
their decay luminosity) trace the centrally peaked den- 
sity profile of the cluster, the average decay luminosity 
would be wiped out by this procedure, and lead to a spu- 
rious depression in measurements of ^. The sky subtrac- 
tion implemented in VIPGI is thus unsuitable for our 
axion search and was not applied. A customized sky- 
subtraction was applied, as discussed in Section flV Bl 

Flux is calibrated separately for each IFU quadrant, 
using observations of a standard star. Finally, the four 
fully reduced exposures are combined. The final data 
cube for A2667 is made of 6806 spatial elements, each 
one containing a spectrum from 3500A to 7000A, and 
covers a sky area of 0.83 arcmin^, centered 5 arcsec south- 
west of the brightest cluster galaxy. The final data cube 
for A2390 is made of 24, 645 spatial elements, each one 
containing a spectrum from 3500A to 7000A, and covers 
a sky area of 3.11 arcmin^, centered 15 arcsec north- 
east of the brightest cluster galaxy. The median spectral 
resolution is ~ 18A. For further discussion of the process 
used to generate the data cubes, see Refs. [47l[50j. 

After producing data cubes in VIPGI, we passed these 
data cubes to a secondary routine that searches for 
emission from axion decay and estimates the noise in 
our spectra. The most obvious source of error is Pois- 
son counting noise. The number of photons observed 
at wavelength A in the j*'' spatial bin is just N\j = 
E[Fx^j/{hc/X)]5tSXSA, where SA = 51.2m^ is the col- 
lecting area of the Melipal telescope, 5X = 5.355A is the 
dispersion of a single VIMOS spectral pixel, St is the 



integration time, F\j is the flux in the j*'' pixel at wave- 
length A, and E is the end to end mean efficiency of 
VIMOS mounted at Melipal. The Poisson counting noise 
is SN\j « ^yN\j, and so SIxj ~ / ^J N\_j. A sec- 
ondary error source is flux contamination from neighbor- 
ing pixels. To include this error, we use the 5% estimate 
of Ref. [43, calculate the 'leakage' contribution to noise 
at each pixel by taking the mean flux of all the nearest 
neighboring pixels, and multiply it by 5%. We also use 
time-logged measurements of the CCD bias and dark- 
current, with the appropriate integration time, to calcu- 
late the additional noise from these sources^ . Finally we 
estimate the flat-fielding noise using the rms difference 
between different sets of efflciency tables. These errors 
are added in quadrature to obtain a data cube of the 
estimated errors in specific intensity. 

In Refs. [H, [l^l, slit locations were chosen to avoid 
the locations of known galaxies, as well as regions that 
showed statistical evidence for faint galaxies 55|- Like- 
wise, we masked out IFU fibers that fell on the locations 
of individual bright sources. Bright sources were identi- 
fied in each cluster image by tagging pixels where the im- 
age intensity exceeded the median by more than la and 
masking IFU fibers that fell on these pixels. Practically, 
this means that 40% of the fibers in each data cube are 
left unmasked. The images used to generate this mask 
are broadband, and so this masking technique will not 
mask out an axion-decay signal. The accidental inclu- 
sion of galaxies could conceivably lead us to erroneously 
attribute their emission to axion decay. This is unlikely, 
given that the spectra of cluster galaxies are dominated 
by continuum emission and line absorption. If we see 
an indication of emission due to line decay, however, we 
may have to revise our masking criteria to take account 
of this possibility. As we shall see later, we imposed up- 
per limits to axion decay, and can safely use the chosen 
masking criterion. The resulting masks were visually in- 
spected to verify that most of the masked fibers fall near 
galaxies. To extract the density dependent component of 
the cluster spectra, we apply a mass map obtained from 
gravitational lensing observations. 



IV. ANALYSIS 
A. Strong Lensing and Cluster Mass Maps 

To model the mass distribution of A2667 and A2390, 
we used both a cluster mass-scale component (represent- 
ing the contribution of the dark-matter halo and the intr- 
acluster medium) and cluster-galaxy mass components as 
in Refs. [i^, [56j. Cluster galaxies were selected accord- 
ing to their redshift (when available, in the inner cluster 
region covered with VIMOS spectroscopy) or their color. 



http: //www . eso . org/observing/df o/quality/VIMOS/toc .html. 
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thus selecting galaxies belonging to the cluster red se- 
quence. For A2667, ISAAC images were used to deter- 
mine J-H colors, whereas for A2390, HST images were 
used to determine I-V colors. The lensing contribution 
from more prominent foreground galaxies was also in- 
cluded, rescaling their lensing properties using the ap- 
propriate redshift. 

All model components were parameterized using a 
smoothly truncat ed p seudo-isothermal mass distribution 
model (PIEMD) [53, which avoids both the unphysical 
central-density singularity and the infinite spatial extent 
of the singular isothermal model. 

The galaxy mass components were chosen to have the 
same position, ellipticity and orientation as their corre- 
sponding images. The K-band luminosity of the galaxies 
was computed, assuming a typical E/SO spectral energy 
distribution (redshifted but uncorrected for evolution of 
constituent stellar populations). Their masses were esti- 
mated using the K-band luminosity, calculated assuming 
a global mass to light ratio [M/L) and the Faber- Jackson 
relation [58*]. The final mass model is made of 70 com- 
ponents, including the large scale cluster halo and the 
individual galaxies. 

Using the LENSTOOL ray-tracing code [111 with the 
HST images, we iteratively implemented the constraints 
from the gravitational lenses. Lensing mass models with 
< 1 were found by fitting the ellipticity, orientation, 
center, and mass parameters (velocity dispersion, core 
radius, and truncation radius) of the cluster scale com- 
ponent, as well as the truncation radius and velocity dis- 
persion of the ensemble of cluster galaxies, using scaling 
relations for early- type galaxies [60| . Cluster galaxy red- 
shifts were measured using the IFU data ^50.) ■ The bright 
central galaxy and several galaxies near the locations of 
strong-lensing arcs were modeled separately from the en- 
semble. The resulting cluster density maps for A2667 
and A2390 are shown in Figs. [3] and 11 [13, Isli . Statisti- 
cal errors in the mass model parameters were propagated 
through the relevant code to produce a fiber by fiber map 
of statistical errors in E. These maps were then used to 
weight different IFU fibers and thus maximize the signal 
to noise ratio of any putative line emission from axion 
decay. 



B. Extraction of One Dimensional Spectra 

Using density maps of A2667 and A2390, we can op- 
timally weight averages over fibers to maximize the con- 
tribution from high density regions of the cluster. This 
maximizes the signal to noise ratio of our axion search 
by emphasizing IFU fibers where maximum signal from 
axion decay is expected. These maps allow us to separate 
emission correlated with the mass profile of the cluster, 
which could be due to axion decay, from a sky back- 
ground that we assume to be spatially homogeneous. Our 
technique is an IFU generalization of the long-slit 'on-off' 
sky-subtraction technique presented in Refs. [l^, [23 |. 
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FIG. 3: Mass map of A2667. The intensity of the image scales 
with density (in units of IQ^^Mq pix~^), where 1 pix = 0.50". 
A density scale is provided on the bottom of the image. The 
horizontal extent of this map is 222.6". The vertical extent 
is 200.0". The thick black line indicates the spatial extent of 
the IFU head on the mass map. 



The real sky background is certainly not perfectly homo- 
geneous, but by making this assumption, we are being 
maximally conservative. With our reduction method, 
any density correlated spatial dependence in the sky 
background will show up as putative emission from ax- 
ion decay. If evidence for a signal is seen, we will have 
to be careful to avoid being confused by sky line emis- 
sion. The projected surface density of the cluster at the 
location in the lens plane associated with the i*'' fiber 
is denoted I]i2,i. Assuming that the only spatially de- 
pendent signal comes from axions, we can then model 
the actual intensity 1^°^ at a given wavelength A and 
spatial pixel i as — (/A/Si2)Si2,i -I- 6a, where 5a 

represents the contribution of a spatially homogeneous 
sky signal, and is the specific intensity in the i^^ 
fiber at wavelength A. Since the signal from axion decay 
is bounded from above by the total component of the 
signal proportional to Si2,i, a measurement of (/a/Si2) 
will either provide evidence of axion decay, or impose an 
upper limit on (/A/Si2)axion- Using a simple linear fit to 
separate the sky background from signal, we extract the 
array (/a/Si2) from each cluster data cube. 

At a small number of wavelengths, this yielded nega- 
tive (unphysical) values for one or both fitted parameters. 
To avoid this, we fit for (/a/Si2) and 5a, subject to the 
obvious constraints (/a/Si2) > and 6a > 0, estimat- 
ing errors a\^i as described in Section UlI CI Estimated 
errors in S12 were also included in the fit. Residuals 
from the best fit are due to fiux noise, imperfect masking 
of galaxies, and variations in fiber efficiency unaccounted 
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FIG. 4: Mass map of A2390. The intensity of the image scales 
with density (in units of 10^^ Mq pix~^), where 1 pix = 0.50". 
A density scale is provided on the bottom of the image. The 
horizontal spatial extent of the map is 157.5". The vertical 
extent is 150.0". The thick black lines indicate IFU pointings 
used to construct our data cubes. 



for by the flat-fielding procedure. As can be seen directly 
from the unconstrained linear-least-squares solutions for 
the parameters (/a/Ei2) and 6a, this procedure places 
higher weights on those VIMOS fibers that fall on higher 
density portions of the cluster. In essence, we used our 
knowledge of the cluster density profile to extract only 
the information that interests us, namely, that part of the 
emission that is correlated with the cluster's mass density 
profile. At some wavelengths, the best fit is (/a/Si2) = 
with very low noise. We verified that these wavelengths 
coincide with those at which a naive linear fit yielded 
negative values for {Ix/T,i2). Thus there is no evidence 
for density correlated emission at these wavelengths. At 
these wavelengths, the emission due to axion decay is 
bounded from above by the brightness of the sky back- 
ground, and so we used (/a/Si2) < bx/i^u) to obtain 
a conservative upper limit on the flux. Decaying axions 
will produce line emission, so it might seem that an ad- 
ditional continuum subtraction might be in order. The 
continuum component of the sky background, however, 
is already subtracted using the techniques discussed, and 
an additional continuum subtraction step would be erro- 
neously aggressive. 

As a test of our sky-subtraction technique, we also 
reimplemented the sky-subtraction technique of Ref. [23| 
and implemented an 'on-off' subtraction by defining 
fibers further than 23" (A2667) or 72" (A2390) from the 
cluster center (defined by the highest density point in the 
density maps) as 'sky' fibers, spatially averaging the flux 
of these sky fibers at each wavelength, and subtracting 
the resulting sky spectrum from each pixel in the 'on' 



cluster region. In this case, sky emission was directly es- 
timated from the data rather than modeled. In this case, 
the best fit for the signal is given by 

where i is a label for the density at the location of a given 
IFU fiber, and a\^i is the error in the specific intensity. 

In the case of A2667, even the fibers furthest from the 
cluster center fall on portions of the cluster where emis- 
sion due to axion decays will be of the same order of 
magnitude as at the center. The sky-subtraction tech- 
nique of Ref. [2^ is thus entirely inappropriate for our 
data on A2667, as it will subtract out a substantial frac- 
tion of any signal and return unjustifiably stringent limits 
to emission from axion decay. 

For A2390, the effective field of view is much larger, 
and so the emission expected from axion decays in the 
outer fibers is much less. Over most of the wavelength 
range of our data for A2390, the different sky-subtraction 
techniques agreed to within a factor of two, leading us 
to believe that our sky-subtraction technique is trust- 
worthy. We used the value for (1^/^12) obtained using 
our sky-subtraction technique, as it is desirable to use 
the same sky-subtraction method for both clusters to 
be self-consistent. Equation ([5]) and the corresponding 
best-fit result in the constrained case essentially yield 
one-dimensional cluster spectra, rescaled by the clus- 
ter density, as shown in Figs. [5] and [H The specific 
intensity values in Fig. [5] were obtained by multiply- 
ing the best-fit values of (/a/Si2) from Eq. ^ by the 

mean (S12) = (E. ^l2,^/<^li) I (E. ^1^1^- The plot- 
ted spectrum is thus not the best-fit spectrum at any 
particular fiber, but an average cluster spectrum. The 
signal to noise ratio of the one-dimensional spectrum ap- 
pears to be higher for A2667 than for A2390, in spite 
of the lower effective fiber number of the data cube for 
A2667. We befieve that this is the case because the data 
cube for A2390 was built using four nights of data, with 
slight variations in sky intensity and efficiency from night 
to night. The subtraction is poorest around the promi- 
nent sky line at 5577A. There is no obvious candidate 
for an emission line due to axion decay. 

C. Limits on the two-photon coupling of axions 

The expected strength of an axion decay line is set by 
Wa,GV through Eq. ([4]), and the expected Gaussian line 
profiles arc shown on top of our appropriately normalized 
upper limits to flux in Fig. [S] for several candidate axion 
masses. The narrow feature at 5 5 77 A, present in both 
panels of Fig. [6l arises from the imperfect subtraction 
of a sky emission line. In the absence of a candidate 
axion decay line, we proceed to put an upper limit on 
the coupling strength ^ of an axion to two photons. 
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FIG. 5: Average one-dimensional sky subtracted spectra 
of clusters A2667 and A2390. Intensity is in units of 
10~^* ergs cm~^ s~^ A ^ arcsec"^. Poorly subtracted sky 
emission lines at 5577A, 5894A, and 6300A have not been 
removed. 
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FIG. 6: Constraints on (/a/Si2) as a function of wavelength 
A for A2667 and A2390. CGS units for specific intensity are 
ergs cm~^ s"""^ A arcsec"^, and E12 — E/(1O^'^M0 pix"'^), 
where E is the projected mass density of the cluster, mea- 
sured using strong lensing. The over-plotted dashed lines are 
theoretical Gaussian spectra for axion decays, with central 
wavelength Aq, corresponding to an axion mass of ma,oV = 
24, 800A(l + z)/Ao. The predicted amplitude is set by Eq. (gl, 
and exceeds the measured values in both the top panel = 
1.0) and the bottom panel = 0.03). 



Since our best-fit values for R\ = {Ix/'^i2) at each 
wavelength come with an error estimate a\, we can cal- 
culate a 95%-confidence limit to the line flux. We assume 
that the distribution of noise peaks is Gaussian, and so 
the probability that an axion decay associated with a 
particular value of i?a,A yields a measured best-fit value 
less than Rx is 



1 



dx. 



(6) 



Eq. ^ yields the 95%-confidence limit on intensity from 
axion decay: 



i?a,A <i?A + 1.65fTA. 



(7) 



At those wavelengths where the best- fit value is i?A = 0, 
we have taken the roughly homogeneous intensity of the 
sky as a very conservative upper limit on the intracluster 
emission. This is many a above the 95%-confidence limit, 
and so at these wavelengths, we just take i?a,A < ^a 
without making our estimate of the upper limit too con- 
servative. Ultimately, we wish to combine the upper 
limits to flux from the two clusters. One of the advan- 
tages of working with two clusters at slightly different 
redshifts is that rest-frame wavelengths falling near sky 
lines (where limits to flux are generally quite poor) at 
one redshift may no longer fall on sky emission lines at 
the redshift of the second cluster. When this is the case, 
we excise wavelengths falling on or near sky lines from 
each spectrum. To account for all the flux in a given 
candidate axion line, in each cluster spectrum, we cal- 
culate the average intensity of non-excised data points 
m a rUa eV window around a series 

of putative line centers spanning the probed axion mass 
range. We weight the noise in the usual way. Assuming 
that our spectra uniformly sample this bin and that flux 
errors are uniform across the bin, we see by integrating 
the Gaussian profile given in Eq. (jH) that 



^a,A — 



2.30 X 10-1* 



'a.oV 



(1 + 2;ci)**5'2(zci)criooo 



cgs. 



(8) 



If axions have the standard thcrmal-freezeout abundance 
[Eq. ([T])], then the limit on the axion coupling is given by 



2.30x10-1** cgs 



1/2 



(9) 



If the cosmological axion abundance takes on some other 
value ilah^, then the limit becomes 



-,1/2 



3.48x10-16 cgs 



(10) 



Since our real bins are not uniformly sampled (because 
of the excision of wavelengths that fall on sky emission 
lines) and since the errors scale with the intensity value 
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at a given wavelength, we make a small correction to this 
expected value. Specifically, 

_ 2.68 X 10-18 emjeY ^jen ^ 

(Aj/(l + z ,)-Aa)^c^ 

where Gj = e ^^^^1 , (11) 

Tx is the set of all non-excised wavelengths lying within 
the bin centered at wavelength A, and j labels wave- 
lengths. The quality of sky background subtraction may 
vary as a result of spatial and temporal variations in the 
sky background from night to night. If it is not due to 
axion decay, the density correlated emission might also 
genuinely vary between clusters. The quantity i?a,A, how- 
ever, will by definition be independent of these factors. A 
simple error weighted mean of the upper limits obtained 
from the two clusters would thus erroneously increase the 
upper limit placed on ^. If two clusters yield different 
best- fit values for Rx, i?a,A must be bounded from above 
by the lesser of these two. By comparing upper limits 
to R\ obtained from A2390 with those obtained from 
A2667 and choosing the lowest value at each wavelength, 
we obtained the maximum values of i?a,A consistent with 
the data. We then applied Eq. (|lip to obtain an upper 
limit on ^ consistent with the spectra of both clusters. 
To account for variation in the upper limits to ^ aris- 
ing from systematic errors in the cluster mass profiles, 
we repeated the preceding analysis, drawing Y^u^i from 
the best-fit NFW (Navarro, Frenk, and White) and King 
profiles to the cluster mass profiles. 

Analytic expressions for the volumetric and surface 
mass density for NFW and King profiles are reviewed 
in Appendix [Xj We determined the mass profile pa- 
rameters (a and a for King profiles, c and a for NFW 
profiles) by fitting to our strong-lensing density maps. 
Using these different density profiles and assuming that 
the mass fraction in axions is = Oa/i^/(r2m^^)j we 
obtain limits to f . The cluster density at the location 
of a given IFU fiber varies from profile to profile, and 
so different fibers receive higher weights when a one- 
dimensional spectrum is extracted. This explains the 
variation in upper limits to ^ that arises when different 
density profiles are assumed. We show the most conser- 
vative (with respect to choice of density profile) limit on 
^ (assuming the thermal-freezeout abundance of axions) 
at each candidate axion mass in Fig.[7l The upper limits 
to ^ in adjacent points along the ma,cV axis are corre- 
lated due to overlapping bins. The narrow black arrows 
near ma,cV = 5.43 and 4.83 mark sharp night-sky lines 
at A = 5577A and A = 6300A, where sky subtraction 
is unreliable and useful limits to ^ cannot be obtained. 
Limits on ^ and specific intensity at the putative line cen- 
ter are displayed for several candidate masses in Table [J 
Our data rule out the canonical DFSZ and KSVZ axion 
models in the mass window 4.5 < ma,oV l£ 7.7, as seen 
in Fig. [51 However, theoretical uncertainties motivate 
the search for axions with values of ^ smaller than those 



TABLE I: Upper limits to central line intensity and ^ at sev- 
eral candidate axion masses, derived directly from sky sub- 
tracted spectra of A2667 and A2390. 



JTla.eV 


(/a/Ei2> {E12) (Cgs) 




4.5 


1.83 X 10"^'' 


7.17 X IQ-^ 


5 


6.04 X 10-^^ 


9.00 X 10-3 


6 


8.74 X 10"^^ 


5.72 X 10-3 


6.5 


9.91 X 10"^^ 


4.60 X 10-3 


7 


9.13 X 10"^^ 


3.41 X 10-3 


7.5 


1.11 X 10"^** 


2.95 X 10-3 


7.65 


8.96 X 10"^^ 


2.47 X 10-3 



allowed by the canonical KSVZ and DFSZ models. 

If we relax the assumption that the cosmological ax- 
ion abundance be given by Eq. ([1]), then our null search 
implies the bound, shown in Fig. [51 to the combination 
^(r^a^^)^^^- We see that if ^ ~ IQ-^, as is the case in 
the KSVZ model, then our results imply an upper limit 
^ah^ ^ 10~* in our mass range, roughly two orders of 
magnitude stronger than CMB/LSS limits (sH], which 
probe densities down to fla,h^ ^ IQ-^. 

D. Revision of past telescope constraints to axions 

As can be seen from Eq. and from the fact 

that the mass fraction of the cluster in axions is Xa = 
^a.h^ / i^mh^), the upper limit on ^ derived from a given 
upper limit on flux depends on the cluster mass model 
used and the cosmological parameters assumed. Refer- 
ences [2^ [131 date to a time when the observationally 
favored cosmology was sCDM (Standard Cold Dark Mat- 
ter: h = 0.5, rim = 1-0, = 0). Moreover, the King 
profiles assumed in those analyses of A2218, A2256, and 
A1413 were based on available x-ray emission profiles of 
the chosen clusters [2^ [l^l (and references therein) . The 
advent of modern x-ray instruments has improved x-ray- 
derived cluster mass profiles, and gravitational- lensing 
studies have allowed measurements of cluster mass pro- 
files, free of the dynamical assumptions required to ob- 
tain a density profile from an x-ray temperature map. 
The quoted upper limits of Rcfs. [H, must thus be 
rescaled, and we have done so up to an ambiguity in slit 
placement for A2218; details are discussed in Appendix [B1 
and the rescaled limits from past work are shown along- 
side our own in Fig . [71 Our limits improve on the rescaled 
limits of Ref. [11, ill by a factor of 2.1 - 7.1. Our fi- 
nal measurement of (/a/Si2) is only noise limited at a 
small fraction (~ 10%) of the available wavelength range. 
The dominant uncertainty is systematic and comes from 
sky subtraction. The expected improvement estimated in 
the introduction assumed that we are limited by Poisson 
noise in the measured fiux, and is thus naive. 

Previous cluster searches for axions used long-slit spec- 
troscopy. Our use of IFU data is novel, and it is conceiv- 
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FIG. 7: Upper limits to the two-photon coupling parameter ^ 
of the axion, derived directly from upper limits to the intra- 
cluster flux of A2667 and A2390. Our data exclude the shaded 
region. The solid and dashed lines show the upper limits re- 
ported in Refs. [23l . [2^ . adjusted (optimistically and pes- 
simistically) for differences between today's best-fit measure- 
ments of the cosmological parameters/cluster mass profiles 
and the assumptions in Refs. [2^, [SJ . Details are discussed 
m Appendix ini The mass range 4.5 eV < ma < 7.7 eV arises 
from the 4000A-6800A usable wavelength range of VIMOS, 
which is smaller than that of the KPNO spectrograph used in 
Ref. [24|. The narrow black arrows near 5.43 eV and 4.83 eV 
mark the sharp night-sky lines at 5577A and 6300A, where 
sky subtraction is unreliable and useful limits to ^ cannot be 
obtained. The shaded exclusion region is derived by applying 
the cluster density profile (strong-lensing map, best-fit NFW 
profile, or best-fit King profile) at each candidate axion mass 
that yields the most conservative upper limit on ^. 
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FIG. 8: Comparison of existing limits to the two-photon cou- 
pling of a 4.5 eV — 14 eV axion with the projected sensitivity 
of our proposed observations of lensing cluster RDCS 1252 
(2 — 1.237). Flux limits and density profiles were assumed 
to be the same as those of A2667/A2390. The best existing 
upper limits to ^ in the higher mass window come from lim- 
its to the Diffuse Extragalactic Background Radiation (DE- 
BRA) , and were rescaled for consistency with today's best-fit 
ACDM parameters and recent measurements [21I . [23 |. The 
limits reported in this and previous work, derived using opti- 
cal spectroscopy of galaxy clusters, are shown for comparison 
24]. Regions inaccessible due to night-sky emission lines 
are marked with narrow black arrows. The two solid horizon- 
tal lines indicate the predictions of the DFSZ and KSVZ axion 
models; the downward arrows indicate that ^ is theoretically 
uncertain. 



able that peculiarities of the data-reduction techniques 
used in IFU spectroscopy may affect the sensitivity of our 
search. To explore this possibility, we have conducted a 
simulation. 



E. Simulation of Analysis Technique 

We simulate axion-decay emission in our data cube for 
A2667, using Eq. ^ and our lensing derived projected 
density maps. We did this at a range of 10 candidate 
axion masses spanning the full mass range of our search. 
We used 3 or 4 different values of ^ at each candidate 
mass. The first value was chosen to be slightly below 
(5—10%) the limit on ^ set by preceding techniques, while 
the second was chosen to be slightly above the upper 
limit. The third and fourth values were chosen to be in 
considerable (factors of 2 and 10, respectively) excess of 
the upper limit. For all simulated axion masses, visual 
inspection of the data cube yields clear evidence for the 



inserted line when ^ exceeds the imposed upper limit. 
An example is shown in Fig. 1101 After inspecting the 
data cubes visually, we applied the routines used for the 
preceding analysis to produce one-dimensional spectra 
for each cube. We then applied the same routine used to 
extract upper limits to ^ to recover the simulated ^ value. 
When the simulated value of ^ exceeded the upper limit, 
we recovered the correct answer in all cases to a precision 
of 5 — 10%. This leads us to believe that our technique 
is robust and our upper limits reliable. References [23, 
[2^ supplement upper limits to ^ derived directly from 
flux limits with limits obtained from a cross-correlation 
analysis. We do the same, using our data on A2667 and 
A2390. 
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FIG. 9; Limits on the combination ^ (f^aft^)^ ^, derived di- 
rectly from upper limits to the intracluster flux of A2667 and 
A2390. Our data exclude the shaded region. Data analy- 
sis proceeds as in the thermal case, but appropriate (more 
general) expressions for the intensity 7a are used. These con- 
straints do not depend on the assumption that axions are 
produced thermally at early times. 



F. Cross-Correlation Analysis 

If there is an emission line at the same wavelength in 
the rest frame of both clusters, the function 

[/J?(x)<ii/;|(i)<ix]"" 

will have a peak at the lag Iq = 
In [(1 -I- ZA2667) / (1 + 2:a239o)], where a; = In A, /i(a;) and 
I2 {x) are the specific intensities of galaxy clusters A2667 
and A2390, and 2:a2667 ^nd 2;a239o are their redshifts 
(23I [2^ . A statistically significant peak in g{l) would 
indicate the existence of an intracluster emission line 
at unknown wavelength (and correspondingly unknown 
axion mass), which could then be searched for more 
carefully in the individual spectra. Peaks due to noise 
may arise either due to the roughly Gaussian fluctuations 
in flux of the individual spectra, or due to imperfectly 
subtracted flux around sharp sky emission lines. It is 
thus appropriate to mask out bright sky lines. If we 
assume that the distribution of remaining noise peaks is 
Gaussian, then the probability that a cross-correlation 
peak with heig ht greater than or equal to s is due to 
noise is [24[ 




1 2 3 4 5 6 7 8 

FIG. 10: The upper left panel of this flgure shows a sim- 
ulated 4255. 2A slice of the A2667 IFU data cube, with an 
axion-decay emission line inserted corresponding to ma,cv = 
7.2 and ^ — 0.011. The flux scale is in units of 
10~^^ erg s~^ cm~^ A ^. This slice, which lies at the ex- 
pected line center, shows evidence of the inserted axion line. 
The resulting 'emission' clearly traces the cluster mass density 
profile. The lower left panel of this figure shows a simulated 
slice of the same data cube, but at 5267. 2A, well away from 
the line center. As expected, no signature of axion emission is 
present this far away in wavelength from the line center. The 
upper/lower right panels of this figure show 4255. 2A/ 5267.2A 
slices, respectively, of the actual A2667 IFU data cube used 
for our analysis. 

Here, tJg is the rms value of the antisymmetric component 
of g{l) and provides an estimate of the correlation due to 
noise, since a Gaussian signal leads to a symmetric cor- 
relation function [2^, [6l[ . Equation fT^]) determines the 
statistical significance of peaks in g{l). Our analysis of 
correlated spectra follows the treatment of Ref. [61] . We 
calculate g{l) using the sky subtracted one-dimensional 
spectra of A2667 and A2390. With a cross-correlation 
technique, we are able to perform a blind search for clus- 
ter rest-frame emission. We find no statistically signifi- 
cant (> 2(Tg) cross-correlation peaks, as shown in Fig.fTT] 
We simulate our cross-correlation based search for an 
intracluster line in order to set alternative upper limits to 
C @, ■ The limits reported in Ref. [23| were obtained 
using this technique. In our simulation, we introduced 
'fake' axion decay lines into both spectra and calculated 
the resulting g{l) for a variety of values ma,cVi thus sim- 
ulating the cross-correlation search for evidence of intra- 
cluster emission. ^ was initially set to a value for which a 
very significant peak in g{l) appeared, and then ramped 
down until the peak at Iq ceased to be statistically signif- 
icant; that is, until P{s) > 5%. At that point, an upper 
limit on ^ was set. The cross-correlation peaks due to the 
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FIG. 11: Cross-correlation function g{l) between sky sub- 
tracted spectra of clusters A2667 and A2390. No peaks with 
the desired statistical significance were seen, at the desired 
lag or elsewhere. A single cross-correlation peak is near the 
expected lag / = 0.00498 for an intracluster emission line com- 
mon to A2667 and A2390. However, it is not statistically 
significant. 



TABLE II: Upper limits to ^ at several candidate axion 
masses, obtained from a simulation of the cross-correlation 
method, using spectra of A2667 and A2390. 



ma,cV 




4.5 


2.73 X 10"^ 


5 


1.08 X 10"^ 


6.0 


9.35 X 10"=* 


6.5 


6.90 X 10"=* 


7 


4.44 X 10"^ 


7.5 


4.31 X 10~^ 


7.65 


2.16 X 10~^ 


axion line arc well fit by Gaussian curves 


throughout the 
1 



rria^cV — (, parameter space. Since the cross-correlation 
function includes the contribution of flux away from the 
line center, rebinning was unnecessary. 

To further distinguish between signal and noise peaks, 
we follow Ref. [l^l in using the usual criterion, |/— < cr, 
where a is the width of the best Gaussian fit to g{l) 
around a cross-correlation peak. We use this criterion in 
both the simulation and the cross-correlation search i23j. 
The resulting limits to ^ at a series of candidate axion 
masses are shown in Table |lTl and are on average a factor 
of '--^ 1.5 less stringent than those derived directly from 
flux. 

One aspect of the correlation analysis of Ref. [23| is 
troubling. Two noisy, imperfectly sky subtracted spec- 
tra were correlated to search for a signal. The analysis 
of Ref. [23j , however, uses one real spectrum (contain- 
ing noise and an imperfectly subtracted sky-background 
signal) with an artificial axion line inserted, and a sec- 
ond, noiseless, template spectrum, containing only the 
artificial axion line, but no imperfectly subtracted sky 
component. Thus, the method simulated in Ref. [2^ 
is not the same as the method used to actually search 
for evidence of an intracluster line, and by artificially re- 
ducing the noise budget of the simulation, could lead to 
artificially stringent constraints. The appropriate way to 
simulate the cross-correlation analysis is to correlate two 
real spectra with artificial axion lines inserted, as we have 
done. Our data also place limits on the decay of other 
relics. 



G. Sterile neutrinos 

Our data might also be used to constrain the decay rate 
of other ~ 5 eV relics, such as sterile neutrinos [4ll. [42. lisl. 
|43 | . Although the prevailing paradigm places the sterile- 
neutrino mass in the keV range, some experimental data 
can be fit by introducing a hierarchy of sterile neutrinos, 
at least one of which is in the 1 — 10 eV range and could 
oscillate to produce photons in our observation window 
62, 63]. In our notation and in the me-.tJ..T ^ "^s limit 
(where is the sterile-neutrino mass), the intensity of 
this signal is 



Sl2 



= 2.4 X 10 



-18 



Bml^y exp - (A, - A,)' cV {2X''y) 



O'l000(l + 2:cl)^<S'2(Zcl) 



■cgs. 



(14) 



where i? is a model-dependent normalization factor, the time abundance of sterile neutrinos [6; 
oscillation is parameterized by a cumulative mixing angle 
9, and the additional power of mass arises from the late- 



flsh^ = 0.3 X TO^^^v sin^ 26*. 



(15) 
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The flux limits in Table |T] impose the constraints 
B < 8.03 X 10-^ 1.14 X 10^'', 9.41 x 10-^ 3.82 x 
10-^ 2.28 X 10"^1.16 X 10"^ 8.12 X 10"^, and 
5.61 X 10^^ for sterile-neutrino masses of nis^eV = 
4.50, 5.00, 5.50, 6.00, 6.50, 7.00, 7.50, and 7.65, re- 
spectively. In conventional models, B — sin'* (2^^)/10^^. 
The parameter B encodes the effects of both the 
early-universe production and the decay of sterile neu- 
trinos, which occurs at the rate Tg^^+'y = 6.8 x 
10-385-1^5^^ sin2 26* [H M, [63, [si. By definition, 
B < 10^*^, and so optical data only constrain sterile 
neutrinos if some novel mechanism increases the oscil- 
lation rate Tg^^^^ by many orders of magnitude. The 
sharp disparity between x-ray and optical constraints re- 
sults from the T cx scaling of the decay rate. 



H. Future Work 



We have demonstrated the utility of applying integral 
field spectroscopy in concert with lensing data to search 
for axions in z ~ 0.2 galaxy clusters. Our technique 
could also be profitably applied to higher redshift galaxy 
clusters. Although flux falls off as I\ (x {1 + Zci)~'^, the 
fact that we are pushing to a higher mass range ma.cv — 
24,800A(H-Zci)/Aa increases the expected signal. Since 
h oc ?i^a,cV' expected signal actually increases as 
I\ (X {1 + Zc\)^ ■ The most distant known lensing cluster 
is RDCS 1252, at redshift z = 1.237 [69]. 

Using existing weak-lensing mass maps for this cluster 
[70| and creating our own strong-lensing maps, we should 
be able to obtain a sky subtracted, spatially weighted 
spectrum of this cluster, attaining flux levels similar to 
those we have obtained for A2667 and A2390. We will 
thus be able to search for emission from decaying axions 
in the mass window 8 eV < ma < 14 eV. Assuming iden- 
tical cluster density proflles and flux limits, we estimate 
the range of ^ values accessible with a telescope search 
for cluster axions in RDCS 1252. The tightest existing 
constraints to decaying relic axions in this mass window 
come from limits to the diffuse extragalactic background 
radiation (DEBRA) [Ml, [13]. As shown in Fig. [3 a VF 
MOS IFU search for axions in this mass window would 
detect very weakly coupled axions, or alternatively, im- 
prove upper limits to ^ by two orders of magnitude. Ap- 
plying Eq. ([1]), we see that 8 eV — 14 eV axions would 
freeze out with abundance 0.12 < < 0.21. An ax- 
ion detection in this mass window could thus account for 
most of the dark matter; a telescope search in this mass 
window would provide a useful check of FSS constraints 
to axion properties. Future discoveries of even higher 
redshift clusters could allow heavier axion mass windows 
to be probed with cluster observations. 



V. CONCLUSIONS 

The axion hypothesis offers attractive solutions to both 
the strong CP and dark-matter problems. A series of 
null searches and astrophysical constraints has narrowed 
down the parameter space of the axion to two mass win- 
dows, one between 10"^ eV and IQ-^ eV, and the other 
between 3 eV and 20 eV. Previous searches for optical 
emission from decaying axions in galaxy clusters have 
constrained the two-photon coupling of the axion in the 
latter window. We have searched for axion-decay light in 
the galaxy clusters A2667 and A2390, taking advantage 
of strong-lensing mass maps of A2667/A2390 to free our 
analysis of dynamical assumptions. Use of the VIMOS 
IFU allowed an increase in effective collecting area, thus 
increasing the sensitivity of our axion search. We ob- 
served no evidence for emission from decaying axions in 
the mass window between 4.5 eV and 7.7 eV. 

Conservatively, we improve on constraints to the two- 
photon coupling ^ of axions by a factor of ~ 3, averaged 
over the entire mass range we explore. This work presents 
the first application of IFU spectroscopy to constrain the 
nature of the dark matter and not just its spatial distri- 
bution. To check that the stringency of our constraints is 
not an artifact of the rather complicated data-reduction 
techniques inherent to IFU spectroscopy, we have simu- 
lated our technique by introducing fake axion lines into 
our data cubes. Our analysis technique accurately recov- 
ers the value of ^, and the axion's signature fades into 
the sky-background as ^ is ramped down below our re- 
ported upper limits. Our simulations demonstrate the 
robustness of our technique, and our work highlights the 
potential of IFU spectroscopy for more sensitive explo- 
ration of the axion mass window between 8 eV and 14 eV. 
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APPENDIX A: KING/NFW PROFILES 

The King profile is parameterized by the expression 

^(^^^ 4.Ga (1-1-^)3/2 ' (^1) 

where a is the cluster velocity dispersion, a is its core 
radius, and r denotes distance from the cluster center. 
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The surface density for a King profile is derived by inte- 
grating by Eq. (|A1[) along the line of sight, and is given 

by 



1 



2^Ga 1 + ^ ' 



(A2) 



where R is the projected radius [71]. The projected mass 
density associated with the NFW mass profile, 



p(r) 



200c^p^Pcrit 



ln(l +CNFW) - J 



cnfw 



+ CNFW 



(A3) 



is E = TsPsfix), where 

fix) = < 




if a; < 1; 

if x = 1; (A4) 



, if a; > 1, 



cnfw is the NFW concentration parameter, and x — 
R/rs ^MMM- 



TABLE B-1: Summary of observations and properties of clus- 
ters used in Refs. 123. [2411. Table entries taken from Ref. [241. 



a 


a 


Inner/Outer aperture 




Cluster (km s~^^ 


[kpc(arcmin)] 


(R/a) 


z 


A1413 1230 


400/1^0^(2.03) 


1.11/4.64 


0.143 






0.65/2.94 




A2218 1300 


200/1^0^(0.88) 


0.94/5.33 


0.171 


A2256 1300 


A73h-o\5.0) 


0.484/2.96 


0.0601 



|;oid ^ _ i]oid(^^g)^ ygi^g ^j^g best-fit values 

assumed in Refs. [1^ [23 . When calculating Snow at the 
slit locations of Refs. [23l [23|. we took the slit locations 
in angular units and obtained physical distances using 
the angular-diameter distance for a ACDM cosmology. 
Applying Eq. jSl, we obtained S^^^^VS^j'^^-i = 0.9853 
and iincw,2/f.oid,2 ^ 449^ ^^^^^ ^^^^ 

seen that this implies ^now,i — l-104^oid,i and ^nGw,2 — 
0.831Coid,2 for A1413. 

The optical depth to lensing by A2256 is very low, 
because of the low redshift of the cluster. As a result, 
lensing derived mass models of this cluster do not exist. 
We took best-fit values from the BcppoSAX x-ray profiles 
of Ref. [77] , in which King profiles are parameterized via 



APPENDIX B: THE EFFECT OF UPDATED 
CLUSTER MASS-PROFILES ON CONSTRAINTS 
OBTAINED FROM A1413, A2256, AND A2218. 

The values of a and a used in Refs. [2^[23] are shown 
in Table IB-1[ along with the relevant redshift values and 
spectral slit locations. In Refs. [H, [13], the sky back- 
ground was removed by subtracting 'off' cluster spectra 
from 'on' cluster spectra. In general, the expected signal 
due to axion decay, in the observer's frame, is 



ho 



Sa(i?)c3 



(Ao/(l + z^j)-Aa)^ 



47rV27rcrAaTa(l + ZclY 



(Bl) 



This can be shown by the same arguments used to derive 
Eq. ([4]). Using this ratio, we can figure out the ratio 
in expected signals. Since I\g cx we can obtain an 
estimate of the upper limit implied by the results of Refs. 
[2^ . [2^ ] , given current measurements of the cluster mass- 
profile and cosmological parameters. 

For A1413, we took best-fit values from the XMM- 
Newton x-ray profiles of Ref. [76(], where it was found 
that A1413 is fit much better by an NFW profile than 
by a King profile. The best-fit NFW parameters are 
Cnfw = 5.82 and r2oo = ?'sCnfw = 1707 kpc. We 



use ^m,ncw^ncw 



0.15, while l^m,oid/ioid — 0-25 is the 



value used in Refs. [23|, |2J]. The projected mass den- 
sity in axions is Sa = [il^h"^ / {^mh^)] S. We define an 
on-off density-contrast S!;*'^ = S!J™(i?on) - S^™(i?off) 
using the best-fit values today. We define another. 



'iR) = 



TcPs 



1 



(B2) 



We then used Eqs. (|B2p and \k2\ to obtain the ratio 
of the best-fit on-off density contrast determined using 
current data to that used in Refs. 123,01: 



sold 



50areCgjp.yy_g 



f^m.now/lL 





1 





























(B3) 



Here, H is the value of the Hubble constant preferred 
today. Using BeppoSax data, best-fit values of cnfw — 
4.57 and Tc = 570 kpc were derived in Ref. [t^ . using a 
redshift of z = 0.0581, and assuming a sCDM cosmology. 
Rescaling this core radius for a ACDM universe, we ob- 
tain Tc = 414 kpc. Inserting these values into Eq. (|B3[) . 
we obtain S!J™/S°id = 0.5982. For A2256, this yields 
Cncw = l-29^oid- If true, recent claims that A2256 is un- 
dergoing merging activity impugn the assumption that 
A2256 is relaxed [zl, [t^]. In that case, the assumption 
of a King profile for A2256 is invalid, and upper limits to 
^ obtained from A2256 have to be revised. 



^ The factor ps usually appears in NFW profiles, and its use in a 
King profile is unusual, but correct. 
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TABLE B-2: Best-fit parameters for the mass model of 
A2218, determined from a strong-lensing analysis. Tlie ta- 
ble was taken from Refs. [s^, [8l[. Square brackets indicate a 
value that was not fit for, but set by hand. The quantity 9 
is the orientation of the ellipse's major axis relative to some 
horizontal in the image plane. 



AR.A.(") 


ADec.(") 


a/b 6{deg) 


rcorc(kpc) 


rcut(kpc) 




+0.2 


+0.5 


1.2 32 


83 


[1000] 


1070 


[+47.0] 


[-49.4] 


1.4 53 


57 


[500] 


580 


[+16.1] 


[-10.4] 


[1.1] [70] 


< 2 


65 


195 


[4.8] 


[-20.9] 


[1.4] [-23] 


< 2 


77 


145 


+0.3 


+0.1 


1.8 53 


< 3 


136 


270 



The strong-lensing analyses of A2218 in Refs. [80|, ISlI 
indicate the presence of several mass clumps in the clus- 
ter, four of which have total masses comparable to the 
total cluster mass, and one, centered on the brightest 
cluster galaxy (BCG), which has a total mass compa- 
rable to a typical galaxy mass. The observed lensing 
configuration is well fit by the set of parameters listed in 
Table [B^ The parameters refer to a PIEMD [57)], whose 
surface mass density is given by 



where a and b are the semi-major and semi-minor axes 
of the best-fitting ellipse, Xc and yc are the best-fitting 
mass centers given in Table IB-2[ translated into physical 
units using the ACDM angular-diameter distance, and uo 
is the velocity dispersion of the cluster. Although these 
lensing data were analyzed using a sCDM cosmology, the 
authors report that the best-fit parameters are insensitive 
at the 10% level to reasonable variations in cosmological 
parameters. The on-off radii are provided without ori- 
entation information in Refs. [l^, [131 j and so we allow 
the slit orientation angle (f> to vary over the full possi- 
ble range, and repeat the preceding analysis to obtain 
a range O.ST^oid < £,now < 0.71^oid- Parameters whose 
values are bounded from above are set to zero for our 
analysis. Depending on the mass bin, the upper limits 
of Refs. [H, [24] come from A2256 or A2218. The up- 
dated upper limits of Refs. [1^ [l^j must thus fall in the 
bracket 0.57foid < ■fnew < l-29^oid- This range is plot- 
ted in Fig. [71 and it is clear that our upper limits to ^ 
improve considerably on those reported in Refs. [23ll2^. 
even when past work is reinterpreted optimistically. 



2G 

^cut 



1 



=2^1/2 





2 

+ 


'y-Vc 


. 1 + e . 




1 - e 



a/b - 1 
a/6+1' 



(B4) 



[1] E. Leader and E. Predazzi, An Introduction to Gauge 
Theories and Modern Particle Physics (Cambridge Uni- 
versity Press, Cambridge, 1996). 

[2] V. Baluni, Phys. Rev. D 19, 2227 (1979). 

[3] R. J. Crewther, P. Di Vecchia, G. Veneziano, and E. Wit- 
ten, Phys. Lett. B88, 123 (1979). 

[4] P. G. Harris et al., Phys. Rev. Lett. 82, 904 (1999). 

[5] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977). 

[6] G. G. Raftelt, Stars as Laboratories for Fundamen- 
tal Physics (The University of Chicago Press, Chicago, 
1996). 

[7] E. W. Kolb and M. S. Turner, The Early Universe 
(Addison-Wesley Publishing Company, Redwood City, 
1990). 

[8] M. S. Turner, Phys. Rev. Lett. 60, 1101 (1988). 
[9] M. S. Turner, Phys. Rev. D 33, 889 (1986). 
[10] D. B. Kaplan, Nucl. Phys. B260, 215 (1985). 
[11] M. Dine, W. Fischler, and M. Srednicki, Phys. Lett. 

B104, 199 (1981). 
[12] A. R. Zhitnitsky, Sov. J. Nucl. Phys. 31, 260 (1980). 
[13] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979). 
[14] A. I. Vainshtein, V. I. Zakharov, and M. A. Shifman, Sov. 
Phys. Usp. 23, 429 (1980). 



[15] D. A. Dicus, E. W. Kolb, V. L. Teplitz, and R. V. Wag- 
oner, Phys. Rev. D 22, 839 (1980). 

[16] G. Raffeh and A. Weiss, Phys. Rev. D 51, 1495 (1995), 
hep-ph/9410205. 

[17] G. G. Raffelt, Phys. Rev. D 41, 1324 (1990). 

[18] G. G. Raffelt, Phys. Rept. 198, 1 (1990). 

[19] G. G. Raflelt and D. S. P. Dearborn, Phys. Rev. D 36, 
2211 (1987). 

[20] K. Zioutas et al., Phys. Rev. Lett. 94, 121301 (2005). 
[21] M. S. Turner, Phys. Rev. Lett. 59, 2489 (1987). 
[22] J. M. Overduin and P. S. Wesson, Astrophys. J. 414, 
449 (1993). 

[23] M. T. Ressell, Phys. Rev. D 44, 3001 (1991). 

[24] M. A. Bershady, M. T. ResseU, and M. S. Turner, Phys. 

Rev. Lett. 66, 1398 (1991). 
[25] E. Zavattini et al. (PVLAS), Phys. Rev. Lett. 96, 110406 

(2006), hep-ex/0507107. 
[26] E. lacopini and E. Zavattini, Phys. Lett. B85, 151 (1979). 
[27] L. Maiani, R. Petronzio, and E. Zavattini, Phys. Lett. 

B175, 359 (1986). 
[28] A. A. Anselm, Phys. Rev. D 37, 2001 (1988). 
[29] M. Gasperini, Phys. Rev. Lett. 59, 396 (1987). 
[30] G. Raffeh and L. Stodolsky, Phys. Rev. D 37, 1237 

(1988). 



16 



R. Rabadan, A. Ringwald, and K. Sigurdson, Phys. Rev. 
Lett. 96, 110407 (2006), hep-ph/0511103. [57 
E. Masso and J. Redondo, JCAP 0509, 015 (2005), hep- 
pli/0504202. [58; 
P. Jain and S. Mandal, Int. J. Mod. Phys. D15, 2095 
(2006), astro-ph/0512155. [59 
T. W. Kephart and T. J. Weiler, Phys. Rev. Lett. 58, [60^ 
171 (1987). 

Y. N. Gnodm, S. N. Dodonov, V. V. Vlasyuk, O. I. Spiri- [61 
donova, and A. V. Shakhverdov, Mon. Not. R. Astron. [62 
Soc. 306, 117 (1999). 

S. Hannestad, A. Mirizzi, and G. Raffelt, JCAP 0507, [63 
002 (2005), hep-ph/0504059. 

G. F. Giudice, E. W. Kolb, and A. Riotto, Phys. Rev. D [64 
64, 023508 (2001), hep-ph/0005123. 
G. F. Giudice, E. W. Kolb, A. Riotto, D. V. Semikoz, [65 
and I. I. Tkachev, Phys. Rev. D 64, 043512 (2001), hep- 
ph/0012317. [66 
M. Kawasaki, K. Kohri, and N. Sugiyama, Phys. Rev. 
Lett. 82, 4168 (1999), astro-ph/9811437. [67 

D. Grin et al., in preparation. 

S. Dodelson and L. M. Widrow, Phys. Rev. Lett. 72, 17 [68 
(1994), hep-ph/9303287. 

A. D. Dolgov and S. H. Hansen, Astropart. Phys. 16, 339 [69 
(2002), hep-ph/0009083. 

X.-D. Shi and G. M. Fuller, Phys. Rev. Lett. 82, 2832 [70 

(1999), astro-ph/9810076. 

T. Asaka, M. Shaposhnikov, and A. Kuscnko, Phys. Lett. [71 
B638, 401 (2006), hep-ph/0602150. 

J.-P. Kneib, R. S. ElUs, I. Smail, W. J. Couch, and [72 
R. M. Sharpies, Astrophys. J. 471, 643 (1996), astro- 
ph/9511015. [73 
O. Le Fcvre et al.. The Messenger 111, 18 (2003). 
A. Zanichelli et al., Publ. Astron. Soc. Pac. 117, 1271 [74 
(2005), astro-ph/0509454. 

T. Moroi and H. Murayama, Phys. Lett. B440, 69 [75 
(1998), hep-ph/9804291. 

S. Tremaine and J. E. Gunn, Phys. Rev. Lett. 42, 407 [76 

(1979). 

G. Covone, J.-P. Kneib, G. Soucail, J. Richard, E. Julio, [77 
and H. EbeUng, Astron. Astrophys. 456, 409 (2006), 
astro-ph/0511332. [78 

E. Julio et al., in preparation. [79 
M. Seodcggio et al., Publ. Astron. Soc. Pac. 117, 1284 
(2005), astro-ph/0409248. [80 
S. D'Odorico et al., The Messenger 113, 26 (2003). 

K. Home, Publ. Astron. Soc. Pac. 98, 609 (1986). [81 

M. Bershady, electronic correspondence. 

G. P. Smith et al., Mon. Not. R. Astron. Soc. 359, 417 



(2005). astro-ph/0403588. 

A. Kassiola and I. Kovner, Astrophys. J. 417, 450 
(1993). 

S. M. Faber and R. E. Jackson, Astrophys. J. 204, 668 
(1976). 

J.-P. Kneib, Ph.D. Thesis (1993). 

P. Natarajan and J.-P. Kneib, Mon. Not. R. Astron. Soc. 

287, 833 (1997), astro-ph/9609008. 

J. Tonry and M. Davis, Astron. J. 84, 1511 (1979). 

A. de Gouvea, Phys. Rev. D 72, 033005 (2005), hep- 

ph/0501039. 

A. de Gouvea and C. Pena^Garay, Phys. Rev. D 64, 
113011 (2001), hep-ph/0107186. 

K. Abazajian, G. M. Fuller, and M. Patel, Phys. Rev. D 

64, 023501 (2001), astro-ph/0101524. 

P. B. Pal and L. Wolfenstein, Phys. Rev. D 25, 766 

(1982). 

A. De Rujula and S. L. Glashow, Phys. Rev. Lett. 45, 

942 (1980). 

K. Abazajian, G. M. Fuller, and W. H. Tucker, Astro- 
phys. J. 562, 593 (2001), astro-ph/0106002. 
V. D. Barger, R. J. N. Phillips, and S. Sarkar, Phys. Lett. 
B352, 365 (1995), hep-ph/9503295. 
P. Rosati, R. della Ceca, C. Norman, and R. Giacconi, 
Astrophys. J. Lett. 492, L21 (1998), astro-ph/9710308. 
M. Lombardi et al., Astrophys. J. 623, 42 (2005), astro- 
ph/0501150. 

J. Binney and S. Tremaine, Galactic dynamics (Prince- 
ton University Press, Princeton, 1987). 
J. F. Navarro, C. S. Frenk, and S. D. M. White, Mon. 
Not. R. Astron. Soc. 275, 720 (1995), astro-ph/9408069. 
J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 462, 563 (1996), astro-ph/9508025. 
J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 490, 493 (1997), astro-ph/9611107. 
C. O. Wright and T. G. Brainerd, Astrophys. J. 534, 34 
(2000), astro-ph/9908213. 

E. Pointecouteau, M. Arnaud, and G. W. Pratt, Astron. 
Astrophys. 435, 1 (2005), astro-ph/0501635. 
S. Ettori, S. De Grandi, and S. Molendi, Astron. Astro- 
phys. 391, 841 (2002), astro-ph/0206120. 
U. G. Briel et al., Astron. Astrophys. 246, LIO (1991). 
M. Sun, S. S. Murray, M. Markevitch, and A. VikhUnin, 
Astrophys. J. 565, 867 (2002), astro-ph/0103103. 
J.-P. Kueib et al., Astrophys. J. 471, 643 (1996), astro- 
ph/9511015. 

G. P. Smith et al., Mon. Not. R. Astron. Soc. 359, 417 
(2005), astro-ph/0403588. 



